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ABSTRACT 

We present the results of a long (~ 93 ksec) XMM-Newton observation of the bright 
BL-Lac object PKS 0548-322 {z = 0.069). Our Reflection Grating Spectrometer (RGS) 
spectrum shows a single absorption feature at an observed wavelength A = 23.33 ± 
0.01 A which we interpret as OVI Ka absorption at z = 0.058, i.e., ^ 3000 kms“^ from 
the background object. The observed equivalent width of the absorption line ^ 30mA, 
coupled with the lack of the corresponding absorption edge in the EPIC pn data, 
implies a column density Aovi ~ 2 x 10^® cm“^ and turbulence with a Doppler velocity 
parameter b > 100kms“^. Within the limitations of our RGS spectrum, no OVII or 
OV Ka absorption are detected. Under the assumption of ionisation equilibrium by 
both collisions and the extragalactic background, this is only marginally consistent 
if the gas temperature is ~ 2.5 x 10^ K, with significantly lower or higher values 
being excluded by our limits on OV or OVII. If confirmed, this would be the first 
X-ray detection of a large amount of intervening warm absorbing gas through OVI 
absorption. The existence of such a high column density absorber, much stronger 
than any previously detected one in OVI, would place stringent constraints on the 
large-scale distribution of baryonic gas in the Universe. 
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1 INTRODUCTION 

Current cosmological models restrict the baryon fraction in 
the Universe to a few per cent of its total matter and energy 
content. A large fraction of that amount of atoms, ions and 
electrons is strongly suspected to reside in the intergalactic 
medium. Lyman a clouds (including Damped Lyman a ab¬ 
sorption systems) are seen to do minate the baryon content 
of th e Universe at high redshift lIStorrie-Lombardi fc Wol^ 
l200(il . but at lower redshifts their number density and the 
subsequent contribution to the ordinary matter content of 
the Universe decrease. 

Detailed simulations of the cosmological evolution of 
the baryons invariably show that the Lyman a absorbing 
gas at temperatures ~ K undergoes shock heating at 
lower redshifts and it s temperature ri s es to 10^ — 10^ K 
ICen^^sUikei ICen et ^ 1200 ll : iFang fc Canizar^ 

DOOoFTDrw^^^i l20(^ . The baryons in the warm and hot 
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intergalactic medium (WHIM) reach 40-60 per cent of the 
total baryon budget in the local universe, according to above 
simulations. Indeed, most of these WHIM baryons live in 
small sparse concentrations, and only those in the deep¬ 
est potential wells (groups and clusters of galaxies) can be 
seen through their X-ray emission. For the remaining of the 
WHIM, resonance absorption lines from highly ionized ele¬ 
ments (OVI, OVII, OVIII, etc.) are the best ch ance to 
detect them iHellsten et allFQQSt iFang fc CanizareslEnOffl . 

Tenuous gas is best detected by absorption rather than 
by emission, if a bright (and featureless) enough source can 
be found, typically an AGN. The limiting equivalent width 
detectable for a resonance absorption line is ultimately de¬ 
termined by the spectral resolution of the spectrograph in 
use (assuming proper channel over-sampling of each resolu¬ 
tion element) but most often limited by the signal to noise 
ratio. Very roughly, for a typical S/N ~ 10 spectrum, where 
each spectral resolution element is sampled by a few chan¬ 
nels, the weakest absorption line detectable has an equiva¬ 
lent width of the order of the width of on e channel. For the 
Refle ction Grating Spectrometer - RGS iden Herder et all 
1200 that we use in this paper, the spectral resolution of 
the first order spectra is around 60 mA, and therefore with 
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a good quality spectrum we might hope to detect lines as 
weak as 10 — 20 mA. With the slightly higher spectral res¬ 
olution grating spectrographs on Chandra, lines as weak as 
5 — 10 mA can be detected, given sufficient signal to noise. 

In the recent years there have been a number of detec¬ 
tions of resonance absorption lines arising in the WHIM, 
both in the soft X-ray band (with Chandra and XMM- 
Newton) and in the Far-Ultraviolet band (with FUSE). 
A local component o f the WHIM was first discovered by 
iNicastro et alJ ^2002^ using Chandra and FUSE d a ta to - 
wards PKS 2155-204. iRasmussen. Kahn fc Paerelsl (l200.'jl 
reported on the discovery of local OVH absorption towards 
PKS 2155- 204, 3 C273 and Mkn 421 with XMM-Newton. 
iFang et*^ ll200.'jl confirmed the detection of a z = 0 OVH 
Kq line towards 3C273 with Chandra. 

The detection of more distant intervening WHI M seen 
in absorption is so far very limited. iFang et all (l2002|l detect 
an OVHI Ka absorber towards PKS 2155-204 at a redshift 
coincident with a n enhanced galaxy densi ty in a well studied 
region of the sky. iMcKernan et alf ll2003^ detect OVHI Ka 
absorption towards the radiogalaxy 3C120, which can be 
interpreted either as intervening gas at z = 0.0147 or as gas 
ejected by 3C120 at a velocity of ~ 5500 km s“^. The current 
status of the detection of feat ures from the non-loca l WHIM 
has been recently compiled bv iNicastro et aJ ll200fJl . In that 
work, it is hinted that WHIM absorption lines can account 
for the missing baryons in the low-redshift Universe. 

In this paper we present ob servations of the BL-La c 
object PKS 0 548-322 z = 0.069 llFosburv fc Disnevll97fll . 
Prior to that iDisnevI lll974^ had detected the presence of a 
group of galaxies at z = 0.04 around this object, suggesting 
that PKS 0548- 322 might be a cluster m e mber . Much more 
recent work by iFalomo. Pesce fc TrevesI (ll995^ found that 
there is indeed a cluster (A S549) around this BL-Lac object, 
but at its average redshift of 2 = 0.069. X-ray emission from 
this cluster is very complicated to detect and measure given 
the brightness of PKS 0548-322. 

We hrst use the RGS spectra to search for absorption 
lines. The single unresolved line found is identified as OVI 
Kq at a redshift intermediate between that of the BL Lac 
and observer. We then use the EPIC pn spectrum to help 
constraining the column density and Doppler velocity pa¬ 
rameter of the absorber. We conclude by examining the tur¬ 
bulence of the absorbing gas and its possible temperature 
under the assumption of ionisation equilibrium. 

2 THE XMM-NEWTON RGS DATA 
2.1 The X-ray data 

PKS 0548-322 was observed by XMM-Newton ll.Iansen et alJ 
I 2 OOIII during its 520th revolution, starting on the 11th of Oc¬ 
tober of 2002, as part of the AO-2 science operations (Obsid 
number is 0142270101). The observation lasted for 94,500 
sec. According to the analysis of the light-curves of the var¬ 
ious X-ray instruments, the radiation environment condi¬ 
tions were good and with low background except for the 
last ~ 10 ks where strong flaring appeared. The pipeline 
products that were delivered to us had been processed with 
version 5.3.3 of the XMM-Newton Science Analysis Software 
(SAS). However, all data products were reprocessed using 


Table 1. Fits to various wavelength regions with a power law 
absorbed by the Galactic column. 


Instrument 

A range (A) 

r 

X^ld.o.f 

^Lots 

RGSl-l-2 

5-35 

1.74 ± 0.02 

2160.26/1351 

- 

RGSl 

5-35 

1.74 ± 0.02 

1051.40/648 

- 

RGSl 

5-15 

2.00 ± 0.12 

196.55/131 

10 

RGSl 

15 - 35 

1.34 ±0.05 

538.07/458 

20 

RGS2 

5-35 

1.74 ±0.02 

1059.15/699 

- 

RGS2 

5-20 

1.98 ±0.04 

599.63/499 

14 

RGS2 

24 - 35 

0.98 ±0.18 

239.10/195 

14 


“ Number of knots used in the spline fit 


version 6.0.0 of the SAS. In addition, the reprocessing incor¬ 
porated a number of effective area corrections with respect 
to the earlier pipeline processed data. 

The RGSl and RGS2 data were filtered for high back¬ 
ground flare intervals, following the data analysis threads 
under the XMM-Newton SOC pages^. This resulted in 
84,600 s of good exposure time for both spectrographs. Both 
spectrographs were operated in High Event Rate mode. 

The spectra of order —1 and —2 of PKS 0548-322 were 
detected in both spectrographs. The 2nd order spectra, how¬ 
ever, had very few counts and were ignored in this work. The 
resulting 1st order spectra (in bins of width around 5 mA in 
the central part of the spectra) have a variable signal to 
noise ratio across the spectrum, which range from virtually 
zero below 6A and above 35A to a peak value of around 
~ 4 — 5 at wavelengths ~ 15 — 20A. 

The search for absorption lines in the RGS spectra, 
needs to start with the determination of the continuum level. 
In order to perform this, we have smoothed the spectra by 
grouping both the RGSl and RGS2 data in bins containing 
at least 50 counts. The regions below 7A and above 35A, 
where the sensitivity of the RGS is very small, have been 
excluded from any further consideration. The determination 
of a good local continuum level is crucial to assess the sig¬ 
nificance of any putative absorption line, and therefore we 
have restricted the continuum fits to patches of the RGS 
spectra where the residuals do not show any large trends. 

To this goal, we have fitted separately the RGS spec¬ 
tra at both sides of the non-functioning GCDs. The RGSl 
continuum has been determined independently in the 7- 12 A 
and 15-35A regions, and the RGS2 continuum has been de¬ 
termined independently in the 7-20A and 24-35A. As it is 
customary, we have used a single power law spectrum ab¬ 
sorbed by a Galactic column of 2.21 x 10^'^cm“^, but in 
reality any smooth function could have been used. (Note 
that in similar optical/UV studies, splines are often used to 
take into account the possible presence of weak broad emis¬ 
sion line co mplexes). The standard fitting routine xspec, 
version 11.2 (lArnaudlll996r was used throughout. The fit¬ 
ted parameters, along with the corresponding values, are 
presented in table Q It is clear that although an absorbed 
power law does not provide a good enough fit to the full 
band, this model delivers a much more acceptable fit when 
restricted to either side of the non-functioning CGDs in both 
spectrometers. We will elaborate further on the continuum 
fitting in Section 2.3. 

^ http://xmm.vilspa.esa.es 
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Figure 1. A portion of the RGSl spectrum, fitted to a power law 
absorbed by the Galaxy, and the ratio of the data to the model. 
The feature at A = 23.33 A is the only significant feature in the 
whole of the RGSl and RGS2 spectra that has a FWHM equal 
or larger than the spectral resolution. 


2.2 Searching for absorption lines 

Detecting absorption lines in a noisy spectrum is a difficult 
task. We have taken the ratio of the data to the contin¬ 
uum model fitted with the binned data, but in the original 
non-grouped format, to search for absorption fines. We have 
chosen the code EQWID, kindly provided to us by Saskia 
Besier at the University of New South Wales (Australia). 
This code searches for significant “connected” regions de¬ 
viating from the continuum, knowing the spectral resolu¬ 
tion. There are only 2 absorption features in the 4 spectral 
patches that have been explored, which are detected with 
significance above 3a. One of them (whose centroid occurs 
at 26.12A) is actually narrower (50 mA) than the spectral 
resolution of the RGSl (60 mA), and is just significant at 
3(7 level. All this suggests that this is a non-real feature 
and in what follows we ignore it. There is, therefore, only 
one absorption line candidate detected by the EQWID soft¬ 
ware, whose centroid occurs at A = 23.33 ± 0.014 A, with 
a FWHM of 74 ± 19 mA and a measured equivalent width 
of 32.3 ± 7.8 mA as returned by that code (1 sigma errors). 
The significance attributed to this detection is ~ 3.7a. Un¬ 
fortunately the RGS2 data does not cover that region, due 
to the failure of the corresponding GCD (number 4) early in 
the mission. 

Fig. El shows a relevant portion of the RGSl spectrum, 
where the absorption fine can be see n. There is a fea t ure in 
the effective area at ~ 23.35 A which iDe Vries et alJ (l200.3ll 
attribute to an ls-2p Oxygen transition in the instrument. 
This feature, which occurs at a slightly, but significantly 
longer wavelength. The wavelength discrepancy between the 
instrumental absorption feature, and the feature we detect 
in our source amounts to ~ 0.02 A, or slightly more than 
one third of an RGS resolution element. The wavelength 
scale has been calibrated to much higher accuracy (8 mA; 
den Herder et al 2001). In addition, as we will show below, 
the measured equivalent width is significantly larger than 
that of the instrumental feature; and the presence of the in¬ 
strumental feature is of course explicitly taken into account 
in the quantitative spectral analysis. We conclude that a 
misidentification of the feature in PKS0548-322 is excluded. 


To further our confidence, we have also inspected other 
high signal to noise RGS spectra of extragalactic objects. 
Rasmussen. Kahn fc Paerelj (l2003ll present RGS spectra of 
3C273, PKS 2155-304 and Mkn 421. In all 3 cases an OVH 
Ka absorption line at z = 0 (21.60 A) is seen, with equiv¬ 
alent width ~ 15 — 26 mA, i.e., weaker than our detection. 
In these spectra, however, there are no features detected at 
around 23.33 A. Section 2.3 presents a further test, based on 
the split of our observations in two parts, which enhances 
the reliability of the detected feature. 

2.3 More on the continuum and significance of 
the absorption line 

All the previous study is based on the assumption that a 
power law absorbed by the Galactic column fits well the 
RGS spectra on either side of the non-functioning CCDs and 
that the calibration (specifically the effective area) is accu¬ 
rate enough. Indeed, the spectrum of PKS 0548-422 could 
have slight, local deviations from the above model and the 
calibration could still be affected by small-scale systematics. 
As it has been mentioned, the way these issues are solved 
in optical/UV absorption line studies is by fitting a local 
continuum -regardless of its physical meaning-, often using 
splines. Visually, the continuum adopted so far around the 
putative feature appears a bit low, and this adds interest to 
this further analysis. 

Therefore we have taken the data to continuum ratio 
resulting from the fits to either side of the non-functioning 
CCDs in both spectrometers and tried to fit a spline to each 
one individually. A region around the absorption fine and the 
instrumental feature ranging from 23.29 A to 23.38 A has 
been masked out from this study. In order to prevent over¬ 
fitting, that could eventually follow statistical noise spikes 
and even potential absorption or emission features, we ex¬ 
plored the decrease of the ^ function of the number of 

knots adopted. We then decided to adopt a number of knots 
below which the not monotonically decreasing. This 

number, an obvious function of the wavelength width of the 
fitted region, is also shown in table Q In all cases the spline 
fitting implies a substantial decrease in the x^ fit to a con¬ 
stant. The data to model ratio for the 4 spectral zones were 
then divided by the spline fits and this was then used for 
further investigation as explained in what follows. 

The routine EQWID was run on the newly re¬ 
normalized ratios. Two features were found with a formal 
significance above 3a. One at A = 14.70 ± 0.02 has a too 
narrow width (FWHM is only 21 mA) to be consistent with 
a real absorption line and it is therefore excluded. The other 
one, is the same fine found without re-fitting the continuum, 
with the same formal significance of 3.7a, the same central 
wavelength A = 23.33 ± 0.012, a similar width of 86 ± 19 
mA and a slightly larger equivalent width of 39 ± 8 mA. 
This latter difference is indeed due to the fact that the con¬ 
tinuum fitted by the splines is slightly larger. However, we 
decided to be conservative and keep the values presented in 
last subsection, bearing in mind that the highly uncertain 
continuum could change the equivalent width of the fine 
probably upwards. 

There is one further experiment we can do with the data 
to model ratio re-normalized to the continua spline fits, given 
the fact that the distribution of these ratios appears approx- 




























































































































































4 X. Barcons et al. 


imately gaussian. This is to assess the real significance of the 
line detected, by means of Monte Carlo simulations. In order 
to perform this, we generated re-normalized data to model 
ratio spectra within the same 4 spectral bands, by assuming 
that values for each channel are gaussian distributed with 
mean 1 and dispersion given by the error bar of the data 
to model re-normalized ratio in the real data. Note that 
since the signal to noise ratio varies significantly across the 
spectra and therefore we cannot use an overall value for the 
dispersion of the gaussian. 

We run 10,000 independent simulations, and on each 
one of them we run the EQWID software. We Hltered out 
those features detected which have either a positive equiva¬ 
lent width (corresponding to an emission feature) or a width 
smaller than the spectrometer’s resolution (60 mA). With 
these Hlters, a total of 605 lines were detected in the 10,000 
simulations with a formal signihcance above 3.7 (t. This in 
practice means that the probability of the detected feature 
at A = 23.33 being a random noise feature is 6%. Indeed 
this implies a real signihcance close to 2a. We have not, 
however, inspected every individual spectrum where a hcti- 
tious absorption line has been detected by EQWID, but a 
random inspection of a few of them shows that many do not 
conform to the shape of an unresolved line and Voigt pro- 
Hle htting either does not converge or results in a very poor 
ht. This implies that the real probability that the detected 
feature is real is actually higher than 2a. 


2.4 Identifying and measuring the absorption line 

To characterize this absorption feature, we have Erst mod¬ 
eled the absorption feature by using the notch model in 
xspec, added to absorbed power law. This notch model de¬ 
scribes a fully saturated absorption “square” feature, which 
is good enough in our case as the line appears unresolved. 
Inclusion of this new component in the ht to the RGSl data 
in the 15-35 A range produces a decrease in the X%.gs from 
1800.17 (for 1620 degrees of freedom) to 1790.38 (for 1618 
degrees of freedom). The F-test yields a signihcance of 98.8% 
to this new component, slightly smaller than the signihcance 
of the line itself, as computed by EQWID. We attribute this 
to the inaccurate modeling of the line by this simple model. 
The best ht is found at a wavelength of A = 23.33!to o 23 A 
and an equivalent width of 32)^^g mA (90 per cent errors for 
one single parameter). Fig. shows the best ht and confi¬ 
dence contours for the observed equivalent width and ab¬ 
sorption line energy using the notch model, both in keV. 

A further test has been performed in order to check 
whether an unhltered bad event (or a set of them) in the re¬ 
gion used to subtract the background could be responsible 
for the absorption feature. Note, however, that the shape 
of the absorption feature conforms to an approximately un¬ 
resolved absorption line. To this end we splitted the RGSl 
event hie in two approximately equal halves in time, ex¬ 
tracted the corresponding spectra, and computed from 
the previously fitted models with and without the notch. 
The residuals of the model without the absorption mod¬ 
eled, both show a negative signal around ~ 23.33 A. Adding 
the notch model, for the hrst half of the observation the 
goes down from 1931.27 to 1926.83 when the absorption 
’’notch” is included with the parameters fitted to the over- 



Figure 2. Best fit (-1- sign) and confidence contours (1, 2 and 3 
sigma) for the line energy and equivalent width of a notch model 
fitted to the absorption lineTriie_staj_sign represents the instru¬ 
mental feature reported by |De^ries_et_ay (2Q23)i which is ~ 2a 
away, and in fact already accounted for in the continuum model. 

all spectrum, and for the second half the y^ also decreases 
from 1895.56 to 1891.64. Therefore, although each half of 
the observation is not sensitive enough for a detection of the 
absorption feature, it is clear that a single bad event in the 
background region is not producing the observed line. 

In order to identify the absorption line, assuming it 
is caused by intervening material, we made an extensive 
search for resonance lines associated to ground configura¬ 
tions, falling in the wavelength range from ~ 23.33 A (in the 
case absorbing material occurs at = 0) down to ~ 21.82 A 
(in the case the line occurs at the redshift of PKS 0548-322 
2 : = 0.069). Note that the most often seen Ka transition 
of OVII (at a rest wavelength of 21.60 A) falls outside this 
range by a large amount. In fact, for the line to be identihed 
as OVII A21.60 the intervening material would have to be 
falling towards PKS 0548-322 at a velocity of ~ 3000 kms“^, 
which is far too large for a cluster of galaxies, especially for a 
poor cluster like A S549. Note that a signihcant velocity dif- 
ference (but in the opp osite sense) has also been found by 
iMcKerua.n et aj J200.3h in an OVIII Ka absorber towards 
3C120. In that case (but not in ours, if the line is OVII 
Kq) there is a possible interpretation of the absorbing gas 
as being in the jet and moving towards us. 

The only physically meaningful option for that line is 
that it is QVI Ka A22.05, wh ich has an oscillator strength 
of 0.576 jPradha.n et al.ir2nn,3h . There is a further line in the 
doublet at A21.87, but its oscillator strength (0.061) is far 
too small to be detectable. Taking that interpretation, and 
the central wavelength measured by the notch model, the 
redshift of the absorbing gas is at z = 0.058 ± 0.001. 

Having established the nature of the line at 23.33 A as 
OVI, it is natural to search for the most often seen OVII 
A21.60 line. In order to do this, we model the RGSl data by 
including two notches tied down at the same (free) redshift, 
but assuming that one is OVI A22.05 and the other one OVII 
A21.60. The ht does not improve with the inclusion of the 
extra notch, but then we have explored a range in parameter 
space to establish limits on the equivalent width of the OVII 
absorption line. Fig. shows the results, hrst in equivalent 
width for OVI and OVII space and then for the OVI/OVII 
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Figure 3. Top: Confidence contours in the equivalent width of 
OVII and OVI parameter space, as obtained by simultaneously 
fitting two notches at the OVI and OVII Kq wavelengths and at 
the same redshift. Bottom: X%,gs with 1, 2 and 3 sigma 

confidence levels) for OVI/OVII and OVI/OV equivalent width 
ratios, as derived from the same data. 


Table 2. Wavelength and oscillator strength of the Ka transi¬ 
tions for the various Oxygen ions considered here 


Ion 

A (A) 

Oscillator strength 

OV 

22.35 

0.565 

OVI 

22.05 

0.576 

OVII 

21.60 

0.694 


equivalent width ratio. The conclusion is that OVI/OVII 
ratio is larger than 4.5, 1.5 and 0.4 at 1, 2 and 3 sigma 
confidence respectively. We have similarly constrained the 
equivalent width of the OV Ka line at 22.35 A. 


2.5 Profile fitting 

Fig. |l]shows the curve of growth for the OVI A22.05 transi¬ 
tion, for a range of Doppler velocity broadening parameters b 
{b = \/2av, (T„ is the velocity dispersion of the gas, assumed 
Maxwellian). Although the range of equivalent widths al¬ 
lowed at 90 per cent level is large, it is clear that some level 



Figure 4. Curve of growth for the OVI transition at 22.05 A. 
The various curves (from bottom to top) correspond to Doppler 
parameters b from 50 to 500 hms”*^ in 50kms“^ steps. The hor¬ 
izontal lines denote the 90 per cent confidence region allowed by 
the data. 


of Doppler broadening needs to be present or otherwise the 
implied value of Aovi would be very large. 

To gain more insight on the physical modeling of the ab¬ 
sorption line, we have performed a Doppler-broadened Voigt 
profile fitting to the ratio of the data to model in RGSl. 
T he Voigt profile has been approximated with the formula 
by IWhitinel (accurate to better than 5 per cent), the 

thermal broadening is parameterized in terms of the b pa¬ 
rameter and the spectrometer response has been modeled 
as a gaussian of FWHM 750kms“^ outside xspec (i.e., the 
model spectrum has been convolved with that gaussian be¬ 
fore being compared to the ratio of data to continuum). Nat¬ 
ural broadening (both radiative and auto-ionising) is only 
relevant to very high column densities and it has been ig¬ 
nored here. Fits were restricted to the 22 — 24 A range. 

Since a fit with the three parameters (Vovi, b and z) left 
free yielded an undefined value for b {b — 435 ± 420 kms“^, 
1-sigma error), we fixed & to a reasonable value of 200 km s“ 
(see later). The best fit provides a substantial improvement 
of the Xvoigt which goes from 182.4 for 171 data points to 
164.94 with 2 degrees of freedom less. The F-test significance 
of this improvement is at the level of 99.98 per cent. The 
best fit values are 2 ; = 0.0580j(g QQQg and log Afovi(cm“^) = 
16.30 ± 0.35 (90 per cent errors). Fig. |21 shows the portion 
of the spectrum fitted, where it is seen that the fit is indeed 
satisfactory. 

There is a previous, shorter XMM-Newton observation 
of PKS 0548-322 (Obsid number 0111830201), performed on 
October 3, 2001. The RGS spectra (which accumulate ~ 48 
ks of good exposure time each ) have been presented and 
analyzed bv iBlustin et alJ (l2004^ which do not find evidence 
for any significant absorption or emission features. Our own 
analysis of that data is entirely consistent with this. In fact, 
including an OVI absorption line with the parameters fixed 
to the above values, to the ratio of data to continuum in the 
22-24A region of the shorter exposure spectrum results in an 
increase of the from 59.13 (for 59 degrees of freedom) to 
66.73. This increase is, however, too modest to exclude the 
possibility that the line is present in the noisier spectrum. 
There is also the possibility of a transient absorption feature. 
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Figure 5. Fit to the ratio of the data to the continuum with 
a Voigt profile for the OVI line, Doppler broadened with b = 
200kms“^. 

as discussed by iBlustin et alJ i2004i for broad absorption 
features in BL Lacs. 


3 THE EPIC DATA 

The two EPIC MOS detectors iTurner et al.l2nnill operated 
in Fast Uncompressed mode, where the central chip (where 
the target was) is read as a single pixel. Since a proper back¬ 
ground subtraction cannot be performed in this mode, we 
ignored the MOS data. 

The EPIC pn detector iStriider et al. I [200111 was op¬ 
erated in Large Window mode. After filtering out high- 
background episodes, the net exposure time where data were 
extracted went down to 77.2 ks. The spectrum of the source 
was extracted, following the SAS analysis threads, using a 
circle of radius ~ 40”. The background spectrum was ex¬ 
tracted from a ~ 5 times larger region free from obvious 
sources, although in a different CCD chip and at a different 
detector y coordinate. 

The count rate from the source region is 15.42 cts“^ 
(without any pattern filtering), slightly above of the count 
rate quoted in the XMM-Newton Users Handbook of 12 
cts“^ where pile-up starts to be relevant. In order to quan¬ 
tify the effects of pile-up, we used the SAS task epatplot 
which extracts the spectral distribution of the events with 
various patterns (singles, doubles, triples and quadruples) 
and compares that to a pile-up free model. The compari¬ 
son shows that in the 0.5-2 keV is very little affected by 
pile-up, with an overall survival of single events of 99.2% 
and an overall survival of double events of 98.7% in that 
band. The differences between the predicted and observed 
spectral distributions do not exceed a few per cent at any 
particular channel in that band, which is of the order of the 
accuracy in the calibration of the EPIC-pn. We then believe 
that the EPIC-pn can be safely used in the 0.5-2 keV band. 
Above that energy, and particularly around 2-3 keV, devi¬ 
ations from the non-piled-up model, exceed 10-15 per cent 
and consequently we ignore these data. 

Thence we have used the 0.5-2 keV band photons, 
keeping only single and double events and only those with 
good spectral quality (FLAG = 0). A fit to a power law 
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Figure 6 . EPIC-pn spectrum and deviations to a power law 
model seen through Galactic absorption. 

with Galactic absorption fixed at Nh = 2.21 x 10^*^ cm“^ 
yields a photon spectral index F = 2.000 zb 0.005 (90 per 
cent errors). The unabsorbed flux in the 0.5-2 keV band is 
S(0.5-2) = (1.865 ±0.004) x 10"“ ergcm-^ s"^ 

Fig- m displays the spectral fit and the residuals, where 
there are no obvious deviations from the model (none ex¬ 
ceeds 5 per cent at any energy in the 0.5 — 2 keV). The fit is 
good, with Xpn = 354.51 for 302 degrees of freedom, leaving 
a probability that the model is wrong of only 2 per cent. Our 
data does not show the previously claimed presence of spec¬ 
tral f eatures in this object as seen by t he Einstein Observa - 
tory llMadeiski et alllioOlTl and ASCA iTashiro et alJlToOlll . 
which were tentatively interpreted as absorption edges. Such 
features, that were once thought to be universal among BL 
Lacs, are absent in the current data. 

We have also tried to search explicitly for the presence of 
OVI in this spectrum, by using the siabs absorption model 
iKinkhabwala et~^ 1200311 . This model keeps track of both 
the photoelectric absorption edge and the resonance absorp¬ 
tion lines (modeled in terms of proper Voigt profiles) for a 
given ion species. Fixing the redshift of the absorber at the 
value measured in the RGS data (z = 0.058) and with a ve¬ 
locity dispersion corresponding to b = 200kms“^ and with 
the RGS best fit Aovi = 10^®’® cm“^ the Xpn improves by 
less than one unit. Fig. Q shows Xpn. a-s a function of Aovi 
for the EPIG-pn data, where it can be seen that the best 
fit values from the RGS data are within the 1 sigma level 
(assuming b ~ 200kms“^). 

4 PROPERTIES OF THE ABSORBING GAS 
4.1 Turbulence 

Since large values of Aovi are not allowed by the EPIC-pn 
data, this poses additional constraints on the curve of growth 
for OVI. In particular, a ~ 30 mA absorption line result¬ 
ing from a large column and very little turbulence (small 
6-parameter) will be inconsistent with the EPIC-pn spec¬ 
trum. 

To quantify this, we have explored a grid of values of 
Vovi from 2 x 10^® to 2 x 10^^ cm“^ and for the Doppler 
parameter 6 ranging from 50 to 1000kms“^. The RGS 1 ratio 
to the continuum has been Voigt profile fitted in the range 
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Figure 7. Xpn (the minimum has been subtracted) as a func¬ 
tion of the column density A^ovi for the fit to the EPIC-pn 
data in the 0.5-2 keV band to a power law absorbed by the 
Galaxy and an OVI absor ber at 2 = 0.058, using the slabs model 
jKinkhabwala et 3|20^. 



Figure 9. Xcomb ^ function of Doppler velocity parameter 5, 
for a joint fit to the EPIC-pn and the RGSl ratio, where the OVI 
column density has been marginalized. The dotted line shows the 
90 per cent confidence level. 



Figure 8 . X^omb ^ function of OVI column density, for a joint 
fit to the EPIC-pn and the RGSl ratio, where the Doppler velocity 
parameter b has been marginalized. The dotted line shows the 90 
per cent confidence level. 

22-24 A and the EPIC-pn spectrum has been fitted in the 
0.5-2 keV with a model consisting of a power law (absorbed 
by the Galactic H column) and a single OVI ion as defined 
by the model slabs at 2 = 0.058. For each value of the 
parameters A^ovi and b the line centre in the RGSl ratio is 
left as a free parameter and the power law in the EPIC-pn 
data is also left as free (2 free parameters). We then minimise 

2 _ 2 I A 2 

Xcomb Xv oigt ' ^Xpn 

Formally the minimum Xcom6 is found at 
log Vovi(cm“^) ~ 16.3 and b ~ 300kms“’^. Fig. |H1 shows 
Xcomb S'® a function of Vovi, coonsidering b an uninteresting 
parameter, from which we derive log Vovi = 16.3 ± 0.3 (90 
per cent errors). Fig.|^shows x^omb ^ function of b, where 
VoVI has been marginalized. At 90 per cent conhdence 
level b = 300^200 (i-6-i b > 100kms“^), but what is 

most important, very small values of b are very unlikely. In 
fact b < 50kms“^ is excluded at > 3 (t confidence, implying 
that the absorbing gas is turbulent. Very large values of b 


cannot be excluded by this analysis because combined with 
smaller values of the column density they can still ht the 
RGS absorption line. 


4.2 Temperature 

To test under what physical conditions can the observed lim¬ 
its on the OVI/0VII and OVI/OV lin e ratios take place, w e 
performed a number of cloudy runs llFerland et a DEigi). 
We assume an optically thin cloud of gas confined by grav¬ 
ity, which is ionised both collisionally and by photoionisation 
from the extragalactic background. We assumed the gas to 
be uniform, with 0.3 solar metallicity, temperature ranging 
from 10"* to 10^ K and particle density ranging from 10“® to 
10“^cm“®. The particle column density was fixed at some 
low but realistic value 10^°cm“^, but line ratios are inde¬ 
pendent on this in the optically thin conditions assumed. 
We adopt the s imple extragalactic parametrization used by 
iNicastro et alJ (l2002ll . as a broken power law with energy 
spectral indices 1.1 and 0.4 respectively below and above 
0.7 keV, norma lized to lOkeV s“^ cm“^ sr“^ keV“^ at 1 keV 
iBarcons et alJl200(il . 

Fig. COl shows the 1,2 and Scr contours allowed in 
the density-temperature parameter space by the measured 
OVI/OVII and OVI/OV line ratios. Note that in display¬ 
ing these contours, the line intensity ratios derived from the 
spectral htting for the various Ka transitions, have been 
converted to column density ratios (which is what cloudy 
computes) by dividing each one by its oscillator strength 
(see Tab. HJ. Collisional ionisation dominates at densities 
> 10“® — 10~'*cm“®, but at the low density end photo¬ 
ionisation by the extragalactic background is the dominat¬ 
ing ionisation process. The 3a upper limit on the gas tem¬ 
perature for OVI being seen but OVII being undetected is 
around 2.5 x 10® K, fairly independent of the density. This 
same value of the temperature is also about the 3 ct lower 
limit for OVI being detected but OV being undetected. 

The conclusion is that it is very difhcult for gas under 
these conditions to show only OVI Ka absorption, without 
OV and/or OVII absorption. 
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Figure 10. Region allowed in the Temperature - Density param¬ 
eter space of the absorbing gas, as obtained from the constraints 
on the OVI/OVII (thin) and OVI/OV (thick) ratios: la (solid), 
2a (dashed) and 3a (dotted). 


5 DISCUSSION 

The XMM-Newton RGSl spectrum of PKS 0548-322 {z = 
0.069) shows a single absorption feature, formally signifi¬ 
cant at the 3.7(7 level. Monte Carlo simulations show that 
this significance is probably lower, but in any case higher 
than 2(7. Intervening OVI Ka absorption at « = 0.058 is the 
most likely interpretation of the absorption line. The equiv¬ 
alent width is ~ 30 — 40 mA, a value that is of the order of 
the one expected to be produced by a group or cluster. We 
do not find any significant absorption feature correspond¬ 
ing to the local WHIM (2 = 0), nor to the putative cluster 
surrounding PKS 0548-322 {z = 0.069). Interpreting the de¬ 
tected absorption line in terms of OVII Ka absorption at the 
cluster redshift, would require the gas inflowing towards the 
BL Lac at a velocity close to 3000 km s“^, which is at odds 
with the expectation that BL Lacs eject material towards 
the observer. 

Despite pr evious claims on the presence of strong ab¬ 
sorption edges llMadeiski et all99li : lTashiro et alll994ri . we 
do not find such edges in the 0.5-2 keV EPIC-pn spectrum 
of PKS 0548-322. Adopting the OVI identity for the absorp¬ 
tion feature and combining the spectral fits to the RGSl 
spectrum normalized to the continuum (Voigt profile fitting) 
with the EPIG pn data, we find that the absorber has an 
OVI column of log Vovi(cm~^) = 16.3T0.3, with significant 
turbulence (Doppler velocity parameter h > lOOkms”^ at 
90 per cent confidence). The lack of OV and OVII Ka ab¬ 
sorption lines in the RGSl spectrum, however, impose al¬ 
most discrepant conditions to the temperature of the ab¬ 
sorbing gas (assumed in ionisation equilibrium by collisions 
and the extragalactic background). Marginal consistency is 
only achieved if the gas has a temperature T ~ 2.5 x 10® K. 

If confirmed this would be the first detection of absorp¬ 
tion by X-ray gas at such a low temperature. Simulations 
do show that absorbing gas blobs in the WHIM s hould span 
a ran ge of temperatures between 10® and 10^ K JPave et all 
I 2 OOIII . although so far only higher ionisation species, proba¬ 
bly corresponding t o the highest temp eratures in that range, 
have been detected. iGhen et aA ll2003r have specifically stud¬ 
ied the cosmological distribution of OVI absorbers via sim¬ 


ulations. Large column density systems, such as the one 
tentatively detected here, are extremely rare according to 
the simulations. If confirmed, this would represent a strong 
challenge to our knowledge of the distribution of baryons in 
Gosmic scales. 

The fact that no other features associated either to the 
intervening absorber itself or to the local WHIM are not 
detected, calls for the need of a better S/N high spectral 
resolution spectrum to confirm the existence of this peculiar 
absorber. 
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